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The heat of combustion of NBS Standard Sample 39i of benzoic acid under standard bomb con- 
ditions has been determined in terms of electrical units. A value of —26,434.0 J g -1 was obtained. 
The total uncertainty in our determination is estimated to be ±3.3 J g" 1 . The uncertainty due to random 
errors was 1.7 J g" 1 and is based on the appropriate factors for the Student t distribution at the 95 
percent confidence limits for eleven determinations of the energy equivalent of the calorimeter and six 
determinations of the heat of combustion of benzoic acid. The principal systematic error, neglect of 
surface temperature correction for our calorimeter, has been assigned a value of ±2.6 J g _1 until more 
reliable estimates of the correction can be made. Particular emphasis was placed on improving the 
precision of a calorimetric measurement over those previously obtained in this laboratory by the use 
of more sensitive auxiliary measuring equipment and more accurate procedures to evaluate the cor- 
rected temperature rise. 
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1. Introduction 

Benzoic acid has served for many years as a ref- 
erence material of known energy of combustion for 
calibrating bomb calorimeters. At the National Bureau 
of Standards (NBS), the Heat Division of the Institute 
for Basic Standards, or its predecessors, has certified 
the energy of combustion of various batches of benzoic 
acid in terms of electrical units for this purpose. 
Maintenance of a capability of a high degree of accuracy 
and precision in bomb calorimetry is a prerequisite 
for this certification activity. 

The calorimetric measurements on which the 
certifications at NBS have been based have been made 
in most instances with a calorimeter of the same 
design, except for slight modifications, as that used 
by Dickinson [l]. 1 In July of 1964, we started a study 
whose goal was to design a new bomb calorimeter 
capable of higher accuracy and precision than existing 
calorimeters. The first step in this project was a study 
of the performance of an existing Dickinson calorimeter 
(NBS S7662). The study was made in the hope that it 
would lead to a better understanding of the limitations 
inherent in the Dickinson calorimeter. 



1 Figures in brackets indicate the literature references at the end of this paper. 



A second reason for making the study is that the 
last certification of benzoic acid in terms of electrical 
units at NBS was made in 1942 [2J. More recent certifi- 
cations have been based upon intercomparison of 
benzoic acid samples. It appeared likely that a cer- 
tification in terms of electrical units could now be 
made with greater accuracy than had been possible 
heretofore, because of the substantial improvement in 
the precision and accuracy of the auxiliary measuring 
instruments of the calorimetric station since 1942. 
Although there was no reason to suspect that the 
intercomparisons of samples had introduced unde- 
tected systematic errors in the certification of batches 
of benzoic acid, we wished to confirm the absence of 
such errors. 

The general features of the calorimeter and its 
method of operation have been adequately described 
elsewhere [3, 4, 5]. A summary of pertinent details 
is given in sections 2 and 4. 

The basis for the correction of the observed tem- 
perature rise of the calorimeter for the effects of heat 
transfer from its environment, Newton's cooling law, 
and of stirring energy is found in the discussion of 
Coops, Jessup, and Van Nes [4]. The particular 
method of calculation used in this work is described 
in section 5 in some detail because of the higher 
accuracy obtained by the use of fewer approximations. 
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Previous analysis of extensive calorimetric data 
from this laboratory had suggested that the uncer- 
tainty in our calorimetric measurements might be 
due mainly to the uncertainties in temperature meas- 
urements during the drift periods. Increasing the 
precision of the measurements of the resistance of 
the platinum thermometer led to a substantial im- 
provement in the precision of measurement of the 
calorimeter temperatures but revealed the presence 
of small systematic deviations from Newton's cooling 
law. Modification of the calorimeter and our method 
of operating it suggested some possible explanations 
for these systematic deviations and reduced but did 
not entirely eliminate them as is discussed in sections 
3, 5.1a, and 6.2. 

A series of electrical calibrations to determine the 
energy equivalent of the calorimeter and a series of 
benzoic acid combustions to determine the energy 
equivalent, in terms of electrical standards, of the 
combustion of benzoic acid are described in sections 
4 and 5. Presentation of the results is completed in 
section 6 by an analysis of random and systematic 
errors. 



2. Experimental Apparatus 

The calorimetric apparatus consists of a stirred- 
water calorimeter surrounded by an isothermal jacket. 
The calorimeter consists of a closed can which con- 
tains a stirrer, a platinum resistance thermometer 
placed near the can wall, a combustion bomb with fuse 
leads and handle, an electrical heater that fits snugly 
around the lower half of the bomb, and a weighed, 
fixed, quantity of water sufficient to be in contact with 
the calorimeter lid after the calorimeter has been 
assembled. The calorimeter temperature is always 
kept below that of the jacket. 

2.1. Calorimeter and Jacket 

The calorimeter and jacket are essentially the same 
as described previously by Jessup [6], except for the 
following modifications. 

The calorimeter jacket has been enclosed by an air 
bath, kept near 27.5 °C with a regulation of ±0.2 °C 
by an on-off controller. The air bath was required to 
make sure the average jacket temperature did not 
change with time. Thermocouple measurements had 
shown the lid (i.e., top) of the jacket was about 0.006 °C 
colder than the rest of the jacket, presumably due to 
poor water circulation in the lid, when the room temper- 
ature was 3 °C below that of the jacket. The jacket 
temperature regulator was replaced by a commercially 
available proportional controller having reset action 
and a nickel resistance thermometer for a sensor. 

To insure a constant stirring rate, the calorimeter 
stirrer was turned by a synchronous motor. 

The oxygen combustion bomb and its internal fittings 
are those described previously [7], except for two 
changes. The fuses were made of 2 cm of 0.002-in-diam 



platinum wire rather than a combustible metal to 
eliminate any energy contribution due to fuse combus- 
tion as described by Prosen [5]. The stem of the bomb 
needle valve was modified so that the inside of the 
bomb could be directly connected to the pressure 
gage of the oxygen manifold during filling, for a more 
accurate pressure measurement. 

The calorimeter heater is an improved version of 
a type described previously [8]. It consists of a 
32-11 heater element of glass-insulated, 0.010-in-diam 
Advance wire which was soft soldered at both ends to 
18-gage, Formvar-covered, copper, current leads. The 
element was inserted in a 3 /i6-in-diam 0.030-in-wall, 
soft copper tube. The tube was flattened against the 
element after the space between the tube and element 
was filled with epoxy-resin cement. The tube (sheath) 
length was selected so that, after the heater was coiled 
to fit the bomb, at least 20 cm of each current lead of 
the heater element was inside the calorimeter. The 
length of each current lead between the calorimeter 
and jacket was 5 cm. Of the two 26-gage copper po- 
tential leads, one was attached to a current lead at the 
calorimeter boundary, the other to the remaining 
current lead at the jacket boundary. 

Four jacket terminals pass through the jacket to 
"temper" thermally the heater leads (or fuse leads in 
the case of a benzoic acid combustion). They are the 
same as those described previously [9], except the 
electrical insulation was changed to a 0.008-in-thick 
layer of Teflon tape and epoxy-resin cement. The insu- 
lation resistance both from the heater element to its 
tube sheath and from the jacket terminals to the jacket 
wall is greater than 100 MH at 100 V. 



2.2. Calorimeter Temperature Measurement 
Equipment 

The platinum resistance thermometer is of the 
Meyers type of construction [10], having a 20-cm-long, 
7-mm-diam Pyrex sheath, and an ice-point resistance 
of about 25.5 ft. The resistance element is wound in a 
single coil very close to the glass sheath to provide 
fast response and to minimize self-heating. For these 
experiments the thermometer head was protected from 
thermal drafts by a cylindrical aluminum shield cov- 
ered with asbestos and aluminum. The thermometer 
leads are connected to the bridge via a selector box. 
The thermometer was reproduceably inserted to a 
depth of 21 cm in the calorimeter can (height, 23 cm). 

The sensitivity of the measurement of the change in 
resistance of the platinum resistance thermometer 
was increased from 10" 5 H (10~ 4 °C) to lO'^ft (10~ 5 °C). 
This was accomplished by replacing the G-2 Mueller 
bridge (smallest dial unit 10~ 4 D) by a G-3 Mueller 
bridge (smallest dial unit 10 -5 ft) and using a more 
sensitive detector of the bridge imbalance. 

The main features of the G-3 Mueller bridge are 
discussed briefly elsewhere [10]. Special shielding of 
the bridge and alteration of the heater supply for the 
bridge thermostat have been described previously [11]. 
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Two galvanometer systems were used with the 
bridge. The new and more sensitive system consists 
of a photoelectric amplifier, having a taunt suspension 
galvanometer, and secondary galvanometer. A V2-mm 
deflection of the secondary galvanometer corresponded 
to a change in thermometer resistance of 1 fifl (1 X 10~ 5 
°C), with bridge current reversal and a thermometer 
current of 2 mA. The less sensitive galvanometer 
system, essentially the same as that described pre- 
viously [12], but with a somewhat more sensitive 
galvanometer, was used to measure the thermometer 
resistance only during the rapid temperature rise of the 
main period of an electrical calibration experiment. 



2.3. Ignition Energy Circuit 

The electrical energy required for ignition of the 
benzoic acid pellet was determined by measuring the 
voltage before (about 25 V) and after (about 24 V) the 
discharge of a capacitor having a measured capaci- 
tance of (40.6 ± 0.5) X 1(F /jlF. The circuit is similar to 
that described by Boyd [13]. 

The fraction of the energy released by the capacitor 
that is dissipated in the part of the circuit external to 
the calorimeter was 0.25 ±0.10. The first step in 
determining this number was to measure the corrected 
temperature rise caused by 50 complete discharges 
of the capacitor, when a short was connected across 
the fuse electrodes in the bomb interior. The fraction 
of energy dissipated external to the calorimeter for this 
particular experiment was calculated using the known 
capacitance, voltages, etc. The fraction dissipated in 
actual combustion experiment was calculated by com- 
bining this data with the measured lead resistances 
and the resistances between the bomb fuse terminals 
in the actual combustion experiments. The main 
source of the estimated uncertainty in the fraction is 
the assumption that the combustion fuse resistance 
does not change before it melts. 



2.4. Electrical Calibration Circuits 

The basic circuit for supplying electrical power to the 
heater has been described elsewhere [4]. 

Two separate, commercially available, Zener- 
diode-stabilized, d-c power supplies were used to 
supply electrical energy to the heater. One unit, with 
a range of 0-2 A and 0-60 V and operated in the 
constant-voltage mode, was used to supply power at 
levels of 60 and 110 W to the heater. The other unit, 
having a range of 0-5 A and 0-105 V and operated 
either in the constant-current or constant-voltage 
mode, was used to supply power at 270 W. The sta- 
bilities of these power supplies in the constant voltage 
mode was 2 to 10 ppm as inferred from measurements 
of heater voltage. 



2 The relative equality and stability (1 yr) of the resistors are ±0.0015 percent and 
±0.0005 percent, respectively. 
•'The other 1,000-il resistor is used only when a volt-box ratio of 1:100 is required. 



The time interval during which power is supplied to 
the heater was measured with a time counter accurate 
to within ±0.0001 s. 

The heater current is determined by measuring the 
voltage across a Reichsanstalt-type, 0.01-11 standard 
resistor. 

The potential drop across the heater is determined 
with a 1 : 1000 nominal voltage divider or volt-box con- 
nected to the heater potential leads. 

The resistive elements of the volt-box are: (a) a 
commercially available unit consisting of ten 2 equal 
10,000-d resistors and two 1, 000-17 resistors hermeti- 
cally sealed in an oil-filled box, and (b) a 100-fl standard 
resistor. Permanent copper links connect the 10,000-fl 
resistors in series with the parallel combination of one 3 
1,000-fl resistor and the standard resistor. Movable 
shorting bars are used to connect the 10,000-(1 re- 
sistors in parallel with this combination. The bars are 
made from copper, have a 1 -cm-square cross section, 
and mercury-wetted surfaces to make electrical 
connections. 

The volt-box is calibrated by measuring the resist- 
ance of the 100-1 ,000-(l parallel combination and then 
the resistance of all the resistors connected in parallel. 
In calculating the volt-box ratio, use is made of the 
fact that, because of the close matching of the 10,000-11 
resistors, the ratio of series to parallel resistance of 
the 10,000-11 resistor set is equal to 100 : 1 within 1 ppm 
[14]. 

The imprecision of the volt-box ratio in 9 determina- 
tions made during the course of the electrical calibra- 
tions was 0.6 ppm (standard deviation of a single 
determination). This is approximately the a priori 
estimated precision of a measurement. The inaccuracy 
was estimated to be between 2 and 10 ppm. 

Potential measurements were made to 0.03 fi\ with 
a six-dial double potentiometer using the more sensi- 
tive of the two galvanometer systems described in 
section 2.2. The range of the potentiometer extends 
to 0.111111 V in steps of 0.1 /jlW. The potentiometer 
voltage reference consisted of three saturated, Weston 
standard cells mounted in an improved version of a 
constant temperature box similar to one described 
previously [15]. Since both short- and long- (3 months) 
term temperature regulation of flie box, as indicated 
by a platinum resistance thermometer in the cell 
compartment, was within ±0.001 °C, the standard 
cell voltage could be assumed to be constant to better 
than 1 ppm during the course of the electrical calibra- 
tions (see sec. 5.3). The potentiometer interdial correc- 
tions were determined in the laboratory before and 
after the set of electrical calibration experiments. The 
potentiometer ratio was checked more frequently 
using an auxiliary circuit similar to that described 
elsewhere [16]. Its constancy, based on eleven determi- 
nations during the course of the electrical calibrations, 
was 1 ppm (standard deviation of a determination). 

The resistors of the auxiliary circuit and the volt-box 
as well as the standard current resistor were kept in a 
stirred-oil bath whose temperature was kept within 
0.01° of 33.16 °C. At this temperature, the rate of 
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change of the resistances with temperature is ex- 
tremely small. Because of the constancy of the bath 
temperature, the resistance of the resistors could be 
assumed to be constant to better than 1 ppm during 
the electrical calibration experiments. 

The resistance ratios required in the calibration of 
the volt-box and potentiometer were measured with the 
Mueller bridge described in section 2.2. 

2.5. Oxygen Handling Equipment 

The oxygen manifold and associated purification 
equipment are essentially as described by Jessup [3]. 
A high-purity grade of commercially available oxygen 4 
was further purified for use in the benzoic acid combus- 
tions. The final oxygen pressure is measured to 0.01 
atm on a calibrated Bourdon pressure gage. The 
temperature of the oxygen in the bomb is measured 
to the nearest 0.1 °C with a calibrated mercury ther- 
mometer inserted in a brass cup that fits snugly around 
the bomb while the bomb is being filled. 

3. Tests of the Apparatus 

Numerous measurements of the variation of the 
thermometer resistance with time during drift periods 
not associated with the principal calorimetric measure- 
ments showed that the imprecision of the measure- 
ments, as indicated by the average deviation of the 
observations from the best smooth curve drawn through 
the data, could be reduced to a few microohms. This 
was done by making resistance measurements at 
1 min intervals with the bridge commutator set in 
alternately the N and then (i.e., next minute) in the 
R positions. Bridge current reversal and interpolation 
to 10~ 6 ft were made in the usual manner [17]. The 
effect of the time variation in the lead resistance was 
eliminated by calculating the resistance R(t), at any 
time t, from the observed readings, adjusted for 
interdial corrections, R' (t) , R'(t-l), R'(t+l) 
according to eq (1). 



R(t) = 



R'(t-D 



R'(t) , R'(t- 



1) 



(l) 



Proper operation of the bridge commutator was 
critical in making these more precise measurements. 
This required replacing the mercury in the switch 
every three to six experiments and periodic reamalga- 
mation of switch contacts. 

Small systematic deviations of the best smooth 
curves drawn through the values of R(t) from the 
resistance-time curves predicted by Newton's cooling 
law were observed. The shape of these deviation curves 
differed. A study of possible causes of the deviations 
was inconclusive. The changes in the calorimeter 
jacket, calorimeter stirrer motor, and the thermal 



shielding of the thermometer head, mentioned in 
section 2, reduced but did not completely eliminate the 
deviations. 

The stirring of the calorimeter was varied in order 
to test the possible effect of variation of the heat 
generated by stirring as a source of deviations. Al- 
though the study of this factor was incomplete, we 
found no clear evidence to indicate that random 
fluctuations of the energy of stirring of the water might 
cause deviations either when the stirrer was operated 
intermittently or when the stirrer was turning at a 
constant rate. The installation of a synchronous motor 
to turn the calorimeter stirrer was made to insure a 
constant rate. We did obtain some results suggesting 
that actual misalinement or flutter in the stirrer shaft 
might be a source of deviations. 

To test the applicability of Newton's cooling law 
over a wider range of temperature than usually occurs 
in a drift period (0.01 °C or less), the average rate of 
change of the resistance of the thermometer was 
measured at six calorimeter temperatures below that 
of the jacket, 12 to 19 resistance measurements were 
made at each calorimeter temperature by the pro- 
cedure given above. Resistances were averaged accord- 
ing to eq (1). The drift rate, AR(t)/At, corresponding 
to a selected thermometer resistance, Rs(t) , was 
determined from the best straight line passing through 
the first order differences plotted as a function of time. 
The number of values of R(t) at each calorimeter 
temperature, the value of Rs(t), AR(t)/ht for the 
selected values of R(t), and the estimated average 
deviation of the first order differences from the straight 
line are given in columns 1, 2, 3, and 4 of table 1. 



TABLE 1. Calorimeter drift rate as a function of calorimeter 
temperature 
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"Obs" 


Av. dev. 




"Obs"-"Calc." a 




10 


28.00717 


568.7 


0.5 






-0.8 




17 


28.06315 


460.7 


0.7 






+ 0.7 




13 


28.11692 


355.0 


0.6 






+ 0.2 




11 


28.16559 


260.0 


0.7 






+ 0.4 




15 


28.21188 


169.0 


0.4 






0.0 




13 


28.25985 


74.7 


0.5 






-0.5 



V~aT/ 



= 1956.35 (28.29828 -/?(*)). 



4 Supplier's impurity analyses were 10.0, 84.0, 0.18, and 15.9 molar ppm for argon, 
litrogen, water, and methane, respectively. 



It may be shown that, if Newton's cooling law holds, 
AR(t)/kt should be a linear function of R(t) with an 
error of less than 3 X 10~ 8 12 min -1 in AR{t)IAt for our 
calorimeter and resistance thermometer. The differ- 
ences between the "observed" values and those cal- 
culated from a least squares fit of the data, assuming 
that R(t) has negligible error in comparison to A/?(f)/Aj, 
are given in column 5 of table 1. The root mean square 
deviation, 0.5 X 10~ 6 O min" 1 , is substantially smaller 
than the value of 3.3 X 10" 6 (1 min 1 obtained pre- 
viously by Jessup [18]. Some of the differences is due 
to the lower sensitivity of Jessup's resistance measure- 
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merit (10~ 5 il) and some is due to the averaging effect 
of eq (l). 5 In any event, the random deviations from 
Newton's cooling law appear to be substantially less 
than those obtained previously. This is consistent with 
an improvement in the calorimetric measurements 
which is discussed further in section 6.1. 

A summary of the magnitudes and types of devia- 
tions from Newton's cooling law during drift periods 
of the calorimetric measurements themselves is given 
in section 5.1. 

4. Experimental Procedures 

4.1. Benzoic Acid Combustions 

The general experimental procedure involved in a 
benzoic acid experiment closely follows that given by 
Jessup [19]. Significant differences are given below. 

The benzoic acid received no special treatment prior 
to the experiment but was used as it came from the 
bottle of NBS standard sample " 39i. Pellets of 
1.51 ±0.015 g were prepared and accurately weighed 
after V2-hr temperature equilibration in a platinum 
crucible in a balance having a readability of 1 /xg. 

After the bomb was filled to a pressure of 30 atm 
of oxygen, 15 to 20 min were allowed to elapse, to 
ensure temperature equilibration of the bomb and 
oxygen, before the final pressure and temperature 
were read. 

The temperature of the calorimeter jacket was 
determined by measuring the resistance of the plati- 
num resistance thermometer to ±10 /jl£1 as originally 
described in section 3 using bridge-current reversal. 

Fifteen to 20 min after the calorimeter had been 
heated to a temperature about 3 °C below that of the 
jacket, measurements of calorimeter temperature were 
started. The preliminary calorimeter-temperature 
measurements with lower sensitivity were essential to 
establish the values of the drift rate and the differences 
between the readings with the commutator in the N and 
R positions as a preliminary step to the much more 
difficult measurements to follow. After 5 to 10 min of 
these readings, the galvanometer sensitivity was in- 
creased to permit measurements to ± 1 /nil, and these 
more precise measurements were carried out for an 
initial period of at least 12 min. 

One minute after the final resistance measurement 
of the initial period, the fuse was ignited by discharging 
the capacitor and initial and final capacitor voltages 
were noted. Time-temperature measurements during 
the rapid temperature rise of the main period were 
made with the thermometer current reduced from 2 mA 
to 0.6 mA. 

After the change in thermometer resistance was 
less than 0.002 (1 min -1 , resistance measurements at 
one minute intervals were jesumed with alternate 
N and R commutator positions at a sensitivity of 
± 10 ixCl. Measurements to ± 1 jjl£1 sensitivity were 



started after the drift decreased to 0.000200 (1 min -1 
and continued through the remainder of the main per- 
iod (defined to end 20 min after ignition) and a final 
period of at least 12 min. The jacket temperature was 
then remeasured as before. 

Table 2 illustrates a typical set of resistance and 
time measurements made during a combustion experi- 
ment. The letters N and R refer to the bridge com- 
mutator settings. Columns labeled average resistance 
are values of R(t) calculated by eq (1) using uncor- 
rected rather than corrected dial readings. The 
course of the experiment with time can be readily 
followed by reference to table 2 while reading the 
foregoing text. 

Table 2. Observations of time and thermometer resistance during a 
combustion experiment 



' Set- section .">. la. 





Resistance minus 28.0 


Average 




Time 
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R 






Min 
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Ohms 


Ohms 


Ohms min ' X 10" 




1 


0.18169 


0.18241 






2 


.18296 








3 




.18368 






4 


.18422 








5 




.18494 






6 


.185487 








7 




.186196 






8 


.186747 




0.186785 




9 




.187451 


.187412 


627 


10 


.188000 




.188037 


625 


11 




.188698 


.188661 


624 


12 


.189246 




.189283 


622 


13 




.189943 


.189906 


623 


14 


.190492 




.190527 


621 


15 




.191183 


.191148 


621 


16 


.191733 




.191769 


621 


17 




.192427 


.192389 


620 


18 


.192967 




.193005 


616 


19 




.193660 


.193622 


617 


20 


.194199 




.194235 


613 


21 




. 194882 


.194849 


614 


22 


.195134 




.195465 


616 


23 




.196110 


.196076 


611 


24 


.196659 




.196686 


609 


25 




.197330 






26 


Ignition 








26.2892 


.20 








26.4802 


.21 








26.5710 


.24 








26.6511 


.26 








26.7207 


.28 








26.7993 


.30 








26.8762 


.32 








26.9740 


.34 








27.0828 


.36 








27.2201 


.38 








27.3731 


.40 








27.4658 


.41 








27.5742 


.42 








27.7227 


.43 








27.8993 


.44 








28.1146 


.45 








28.2697 


.455 








28.4495 


.460 








28.6790 


.465 








29.0040 


.470 








29.1801 


.472 








29.3914 


.474 








29.6697 


.476 








29.8461 


.477 








30.0682 


.478 








30.2112 


.4785 








30.3543 


.4790 








30.4885 


.4795 








30.7058 


.4800 








30.9726 


.4805 








31.1630 


.4808 








31.4654 


.4812 








32 


.48175 








33 




0.48230 






34 


.48241 








35 




.482614 
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Table 2. 



Observations of time and thermometer resistance during 
a combustion experiment — Continued 



Time 


Resistance 


tiinus 28.0 


Average 


Drift 






resistance 




N 


R 


minus 28.0 




Min 


Ohms 


Ohms 


Ohms 


Ohms min~ l x 10« 




- 








36 


.48258* 








37 


.482684 








38 






39 




.482830 






40 


.482769 








41 




.482915 


0.482863 




42 


.482854 




.482906 


43 


43 




.483000 


.482948 


42 


44 


.482938 




.482989 


41 


45 




.483080 


.483031 


42 


46 


.483024 




.483074 


43 


47 




.483166 


.483116 


42 


48 


.483107 




.483158 


42 


49 




.483250 


.483200 


42 


50 


.483192 




.483241 


41 


51 




.483331 


.483282 


41 


52 


.483274 




.483323 


41 


53 




.483413 


.483364 


41 


54 


.483355 




.483405 


41 


55 




.483496 


.483446 


41 


56 


.483437 




.483488 


42 


57 




.483580 


.483530 


42 


58 


.483522 




.483572 


42 


59 




.483664 


.483613 


41 


60 


.483603 




.483654 


41 


61 




.483746 


.483695 


41 


62 


.483684 




.483736 


41 


63 




.483830 


.483778 


42 


64 


.483769 




.483819 


41 


65 




.483908 


.483859 


40 


66 


.483852 




.483900 


41 


67 




.483989 


.483941 


41 


68 


.483934 




.483983 


42 


69 




.484076 


.484025 


42 


70 


.484013 









In this way it was possible to obtain up to eight meas- 
urements of both the heater current and voltage during 
a 5 min heating period. 

At frequent intervals during the course of the calibra- 
tion experiments, the volt-box was calibrated, the 
potentiometer ratio was determined, and the constancy 
of the resistance of the 0.01-d standard resistor was 
checked. 

5. Calculation of Results of Calorimetric 
Experiments 

A series of 6 benzoic acid combustion experiments 
and then a series of 16 electrical calibration experi- 
ments were carried out. The electrical calibration 
experiments are separated into three groups depending 
upon whether the amount of power supplied to the 
calorimeter heater was approximately 110, 60, or 
270 W. 



5.1. Calculations of the Corrected Temperature Rise 

Values of corrected temperature rises are listed for 
the appropriate experiments in tables 3 and 4. Ob- 
served temperature rises were of the order of 2.82 °C, 
and the correction due to heat transfer between the 
calorimeter and jacket was about 1.5 percent of this 
value. 



4.2. Electrical Energy Measurements 

Much of the experimental procedure was similar to 
that outlined for the chemical heat measurements in 
the previous section. The procedure for making the 
electrical calibration measurements is essentially that 
outlined previously [4, 2]. 

The combustion bomb was prepared without a pellet 
or platinum fuse but contained the usual platinum 
crucible and the usual amounts of oxygen and water. 

Just prior to a run, the temperature of the standard 
cell enclosure was measured. The emfs at the output 
leads of the volt-box and the 0. 01-11 standard current 
resistor were measured before and after a run when no 
power was supplied to the heater to obtain values of the 
residual (thermal) emf of the circuit. 

Thermometer-resistance measurements during the 
rapid temperature rise (i.e., while the electrical power 
was on) of the main period were made in the same 
manner as in the benzoic acid experiments. Measure- 
ments of the potentials at the output of the volt-box 
and across the current resistor which were 0.054 and 
0.018 V, respectively, for the 110 W calibration experi- 
ments were made between temperature measurements. 
(The time of each potential measurement was differ- 
entiated from those of the main period temperatures 
by actuating the timer-printer with the standard second 
signal to record two times, exactly 1 s apart, which 
bracketed the time of the corresponding measurement.) 



a. Treatment of Initial and Final Period Data 

The observed readings of the bridge during the 
initial and final periods, illustrated in table 2, were 
converted to ohms by applying the various bridge 
corrections and then averaging according to eq (1). 
Temperatures corresponding to each R(t) were calcu- 
lated by the Callendar equation [10]. The parameters 
in this equation that are characteristic of our thermom- 
eter were determined by the Temperature Section of 
the Heat Division, NBS. 

The calorimeter temperatures, 6(1, J) , for the initial 
{1=0) and the final (1=1) drift periods were fitted to 
the integral form of Newton's cooling law, eqs (2). 



d(I,J)=C(I)+D(I)-P(J) 
P(J) = (l-e- KJ )K. 



(2a) 
(2b) 



In eq (2), C(I) and D(I) are different constants for 
each drift period, J is the time of occurrence of the 
temperature relative to the start of the appropriate 
drift period, and K, the cooling constant of the calorim- 
eter, is defined by eq (3). 



K=[D(0)-D(l)]/[C(l)-C(0)]. 



(3) 



A consistent fit of the temperature-time data by 
eqs (2) and (3) was started by estimating values of 
C(I) ,/)(/), and K from the first and last datum points 
of each drift period. Calculation of "improved" values 
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Table 3. Benzoic acid combustion calculations 



Expt. 
No. 


Mass 


FxlO K 


M 


nhd+FVM 


E' (cont) 


A£k;n 


-A£(HN0 3 ) 


A* 


E (cal) 


















(-A£ /( ) 




g 




°C 


g*C ' 


y°<: ' 


J 


J 


(-g°C-')xlO K 


«r°C-' 


1 


1.523394 


- 189.8 


2.849658 


0.5344869 


+ 1.43 


+ 1.49 


+ 0.49 


-27.8 


0.5344591 


2 


1.507283 


-166.7 


2.819376 


0.5345267 


+ 1.79 


+ 1.60 


+ 0.52 


-39.3 


0.5344874 


3 


1.492242 


- 166.9 


2.791449 


0.5344869 


+ 1.80 


+ 0.93 


+ 0.70 


-46.0 


0.5344409 


4 


1.493978 


- 160.8 


2.794708 


0.5344880 


+ 1.88 


+ 1.23 


+ 0.58 


-46.6 


0.5344414 


5 


1.525560 


-157.2 


2.853627 


0.5345198 


+ 1.93 


+ 1.21 


+ 0.58 


-46.3 


0.5344705 


6 


1.501083 


- 180.8 


2.807857 


0.5345043 


+ 1.57 


+ 1.38 


+ 0.52 


-33.8 


0.5344705 



Average 

Std. dev. of experiment.. 
Std. dev. of mean 



0.5344616 

0.0000183(0.0034%) 

0.0000075(0.0014%) 



rA£ I( ; N -A£(HNO :t ) 
(+A0) 



E 1 (cont) 



(-A£«) 



Table 4. Calculation of the energy equivalent of the calorimeter 



Kxpt. 


1 EvE,dt 


MFi) f E*dt 


QrlFt 


Qe 


M 


— E' (cont) 


E (cal) 


No. 


( P 2 s)xl0--' 


(p 1 8) x 10 ,:t 


(p 2 s) x 10 +: < 


J 


°c 


y°c-« 


J °c ' 









Heater po 


*er 110 W 








7 


3550.0634 


1.1258 


0.0377 


38.955.397 


2.757342 


0.30 


14,128.18 


8 


3559.4812 


1.1288 


.0377 


39.058.739 


2.764372 


0.32 


14,129.66 


9 


3615.7432 


1.1467 


.0378 


39,676.104 


2.808623 


0.39 


14,126.92 


10 


3605.5325 


1.1434 


.0378 


39,564.063 


2,800523 


0.02 


14,127.40 


11 


3585.7557 


1.1372 


.0378 


39,347.050 


2,785424 


0.20 


14,126.25 


12 


3591.6315 


1.1390 


.0378 


39,411.526 


2.789739 


0.16 


14,127.48 


13 


3578.4748 


1.1348 


.0378 


39,267.158 


2.779515 


0.16 


14,127.50 


14 


3591.6153 


1 . 1 390 


.0378 


39,411.348 


2.789528 


0.14 


14,128.46 


15 


3588.6870 


1.1381 


.0378 


39,379.216 


2.787293 


0.14 


14,128.26 



Heater power 60 W 



3575.0618 
3573.4132 



1.1342 
1.1337 



0.0069 
.0070 



39,229.362 
39.211.273 



2.776922 

2.775357 



0.07 
0.07 



14.126.99 
14,128.44 



Heater power 270 W 



18* 


3582.3213 


1.1348 


7.3177 


39,389.320 


2.788869 


0.13 


14,123.89 


19* 


3548.3917 


1.1240 


7.3241 


39,017.047 


2.762050 


0.13 


14,126.25 


20* 


3506.1535 


1.1107 


7.2380 


38,552.619 


2.729664 


0.16 


14,123.74 


21* 


3588.2289 


1.1368 


7.2382 


39,453.237 


2.792813 


0.16 


14,126.86 


22 


3612.1828 


1.1440 


0.0208 


39,636.866 


2.805518 


0.19 


14,128.36 



nt mode of power supply; all others eonstant voltage modi-. 



of these constants proceeded by computing the devia- 
tions, F(7, 7), of the "observed" temperatures from 
those calculated by eq (2) using the initial estimates of 
C(/), D(I), and K. These deviations were fitted by 
least squares with deviation equations linear in 
P(J) according to eq (4). 



F(/,y) = B(o,/)+fl(i,/)P(./). 



(4) 



The values of fl(0, /) and B(l, /) were used to correct 
the initial estimates of C(/), D(I) , and K. The cycle of 
calculations, performed by a high speed digital com- 
puter, was repeated until either the sum of the squares 
of the deviations of the initial and final drift periods 
increased or the values of K calculated on two succes- 
sive cycles differed less 6 than 1 X 10 7 min" 1 . 



•The latter "stop" was inserted because there is no assurance of a distinct minimum in 
the sum of the squares of the deviations. It was. in fact, the terminating step of all the cal- 
culations alter two or three iterations. The value of 1 X 10 7 min ' is based on the fact that 
the jacket temperature was constant to only about 1 X 10 4 °('., and A. is approximately 
2X 10'' min '. 



The observed temperature rise was computed from 
the calculated values of initial and final temperatures 
using eqs (2) and the final values of C(7), D(I), and K. 

This method of curve fitting assumes the values of 
6(1, J) of each drift period are independent which 
cannot be so since the temperatures are calculated 
from observed resistances which were averaged ac- 
cording to eq (1). 

It can be shown that, because of the use of eq (1), 
runs or deviations of the same sign should occur when 
the difference between "observed" and calculated 
temperatures are plotted as a function of time even 
though the errors are random. Usually, runs of three 
to five deviations of the same sign were observed. 
This indicates that in the case of random errors, those 
in the resistance measurements predominate. The 
systematic errors in the values of "observed" tem- 
perature and the calculated temperatures caused by 
using eq (1) are negligible. 

Although the deviation plots sometimes showed 
systematic deviations from Newton's cooling law, they 
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were small. The average root mean square deviation 
of "observed" from calculated values of 6(1, J) was 
slightly less than 2 X 10" 5 °C. 



b. Treatment of Main Period Data 

The correction to be added to the observed tem- 
perature rise due to heat transfer from the jacket and 
stirring energy is calculated according to a method 
discussed by Dickinson (see [20, 4]). In this method, 
the main period data is used to evaluate the parameter 
t x , defined by eq (5), and the remaining parameters 
are determined from the values of C(7), D{I), and K 
determined from the drift period data in section 5.1a. 



tx=t { 



(Of 



Si) Jr. 



(6-6i)dt. 



(5) 



In eq (5), and 6j are the temperatures at the begin- 
ning, h, and end, t/, of the main period, respectively. 

To keep the error in the calculation of the correc- 
tion to the observed temperature rise to 1 X 10 -5 °C 
in our experiments, t. r must be calculated to 0.1 s and 
the integral in eq (5) must be calculated with an ac- 
curacy of 1 part in 10,000. The following two procedures 
were adopted to minimize the calculation error. 

First, the temperature-time data were graphed to 
eliminate gross errors and fill in gaps where data was 
sparse by superimposing corresponding portions of 
the main period of similar experiments. 

Second, because eq (5) may be rewritten with an 
integral whose integrand is independent of tempera- 
ture, temperatures of the main period were calculated 
in the same manner. Since virtually all the resistance 
measurements during the rapid temperature rise 
were made with fixed commutator position, the formula 
for the estimated correction for thermometer lead 
resistance applied to these readings was used through- 
out the main period. Thus, although the method of 
calculation of temperatures was in other respects the 
same as that for the initial and final periods, 6; and 6/ 
in eq (5) are not quite the same as those used in cal- 
culating the observed temperature rise. We assumed 
the thermometer lead resistance correction varied 
linearly with the thermometer resistance in the main 
period from the average initial to the average final 
drift correction. Since the difference between these 
corrections was about 15 fxCl for all the experiments, 
the net effect of an error in the assumption is negligible 
in comparison to other errors in the treatment of the 
data of the main period. 

The integral in eq (5) was evaluated using the spline 
numerical integration procedure whose character- 
istics are described in detail elsewhere [21]. Essen- 
tially, the method fits a piecewise cubic with continuous 
first and second derivatives to the data. Calculations 
were carried out with a high speed computer using a 
program written in the Dartmouth BASIC language 
[22 j. This program was based on a subroutine origi- 
nally written in FORTRAN [23]. 



5.2 Benzoic Acid Combustion Calculations 

Since benzoic acid is certified as a secondary heat 
standard for use in determining the energy equivalent 
of bomb calorimeters, it is convenient to express its 
energy of combustion per unit mass of benzoic acid 
under so-called standard bomb-conditions [4]. Deter- 
mination of the energy of combustion of benzoic acid 
under standard bomb conditions, Mb, is carried out 
in two steps, the first of which is summarized in this 
section. 

In order to compare the benzoic acid combustion 
data with the electrical calibration data, we defined the 
standard calorimeter to be the calorimeter can con- 
taining sufficient water so that their combined weight 
is 3750 g, the calorimeter lid, the electrical heater, 
and the combustion bomb with its standard contents, 
its external fuse leads and its handle. The standard 
bomb contents were defined to be 1 g of water, the 
platinum crucible, platinum fuse, and 0.441 moles 
of oxygen (30 atm at 25 °C). 

By selecting the temperature to which the iso- 
thermal bomb process is referred to be the final tem- 
perature the calorimeter would have if the heat transfer 
from the calorimeter environment and stirring energy 
were zero, the ratio of the energy equivalent of the 
standard calorimeter, ^(cal), in joules per degree, to 
M B , in joules per gram, can be calculated from the 
experimental results by eq (6) (see [20], [5]). 

£(cal) (l + F)m s 



M B 



+ 



A0 



M lGN -M h 



(A0) 



-E' (cont) 



1 



J (-M B ) 



(6) 



M b '=(1+F)Mb 



(7) 



In eq (6), Ad is the corrected temperature rise, m s is 
the mass (i.e., weight in vacuo) of benzoic acid in 
grams, M Kj ^ is the electrical energy supplied to the 
calorimeter to ignite the benzoic acid, M Hm:i is the 
energy of formation of nitric acid produced in the actual 
combustion reaction, and E' (cont) is the correction 
added to E(ca\) to give the energy equivalent of the 
actual calorimeter prior to sample combustion. The 
quantity F is related to the energy of combustion, 
M' H , per gram of benzoic acid under actual bomb 
conditions by eq (7). 



Numerical values of F were calculated from the ap- 
proximate formula given by Coops, Jessup, and Van 
Nes [4]. 

Since the second term on the right-hand side of eq 
(6) constitutes only a small contribution to the total, 
an approximate value of M H may be used in this term. 
For this we used —26,434 J g 1 . An error of 10 J g _1 in 
this value of M B would produce an error of only 1 ppm 
in-£(cal)/A£„. 
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A summary of the calculations is given in table 3. 
In calculating F (col. 3), the pressure of oxygen at Oi, 
was calculated from observed pressures using the 
equation of state for oxygen given by Coops, Jessup, 
and Van Nes [4^. £'(cont) ? _ listed in column 5, was 
evaluated using 1.21 J g _1 °C -1 [4] for the heat capac- 
ity of benzoic acid and a value of 21.17 J mol -1 [20] for 
the constant volume heat capacity of oxygen at 30 atm. 
Measurement of AEicx is summarized in section 2.3 
and the determination of A£'(HNO:0 was made in the 
usual manner [4] . 

More accurate values of F were calculated by eval- 
uation of the Washburn correction by the procedure 
described by Hubbard, Scott, and Waddington [20] 
using the computer program developed by Shomate 
[24]. A very slight difference was observed in the F 
values, for which a correction is made in the calcu- 
lation of (—AEb) at the bottom of table 5. (The dis- 
persion of the values of E (cal) I (— AE h) was not 
altered significantly.) 

5.3. Electrical Calibration Calculations 

The energy equivalent of the calorimeter, £(cal), 
is related to the total electrical energy, Q e , supplied 
to the calorimeter during a calibration experiment, by 
eq(8) 



Q e = [£(cal)+F(cont)]A0. 



(8) 



In eq (8), AS is the corresponding corrected tempera- 
ture rise and E'(cont) was defined in section 5. Q<> 
was evaluated by eq (9) 



Qe = F x P E v E,dt-F 2 I ' Efdt + Qr 



(9) 



In eq (9) the various terms on the right are, respec- 
tively, the energy calculated from the potentiometer 
measurements, the correction for current passing 
through the volt-box, and correction for transients 
(i.e., any voltages not measured with the potentiom- 
eter). Ey and Ej are the observed dial readings of the 
potentiometer at the output of the volt-box and across 
the current resistor, respectively after applying inter- 
dial and "zero" corrections. T is the length of time 
power is supplied to the heater. F\ and F> are constants 
which depend upon the volt-box ratio, resistance of 
the current resistor, the standard cell voltage, and the 
calibration of the potentiometer. 

The numerical values of the terms on the right of 
eq (9) divided by F\ are given in the first three columns 
of table 4. The unit p in these tables stands for a 
"potentiometer unit" and is nominally 1 V. The values 
of Q ( >, A0, £'(cont), and £(cal) are listed in the remain- 
ing columns of table 4. 

Because of the constancy of the volt-box and poten- 
tiometer ratios and the inferred constancy of the 
standard current resistor and standard cell voltage to 
1 ppm or better, F\ and F> were assumed to be con- 
stants. This assumption causes negligible imprecision 
in £(cal) as compared to that from other sources. The 



numerical values of F\ and F> are 1.09765196 VKl~ l pP 
and 1.099975 X 10~ 4 , respectively. 

The standard cell voltage was determined by the 
NBS Electrochemistry Section of the Electricity Divi- 
sion by comparison with the NBS voltage standard. 
A correction of 1.4 microvolts [25] was added to this 
calibrated value since the average temperature of the 
cells during our measurements was 0.024 °C less than 
at the time the cells were calibrated. 

The resistance of the standard current resistor was 
based on its certified resistance at 25 °C and the tem- 
perature coefficient of its resistance as determined by 
the supplier. The maximum error in this calculated 
value of resistance at the oil bath temperature (33.16 
°C), as determined by ten potentiometric comparisons 
during the electrical calibration experiments of its 
resistance with that of a standard 0.1 (1 resistor cali- 
brated by the NBS Resistance and Reactance Section, 
was found to be 0.002 percent or less. After our meas- 
urements were completed, a more accurate value for 
the resistance was determined by the NBS Resistance 
and Reactance Section. The corresponding correction 
to the data was made in the final calculation of (— AEb) 
in table 5. 

The first potentiometer measurements were made 
some 10 to 25 s after power was applied to the calo- 
rimeter heater. A total of 6, 13, and 2 measurements 
of Ev or Ei were made for the 110,60, and 270 W experi- 
ments, respectively. Smoothed values of E\ and E, 
were obtained from separate plots of the observed 
values as a function of time; the time dependence, 
more or less unknown, of the 270 W measurements 
had to be inferred from other considerations [26]. 
Values of the integrals in eq (9) were computed with 
the spline integration subroutine of the computer 
program used in the calculation of the integral in eq 
(5). 

After the electrical calibrations were completed it 
was found that large voltage transients occur when the 
power source supplying 270 W to the calorimeter 
heater in the constant current mode is switched from 
the dummy heater to the calibration circuit. These 
transients had disappeared well before the first meas- 
urements of E\ and Ei were made with the poten- 
tiometer. Subsequent examination of the voltage-time 
characteristics of this power source in the constant 
voltage mode and of the power source used to supply 
110 and 60 W to the calorimeter heater in the constant 
voltage mode also indicated the presence of tran- 
sients, but these were much smaller. 

Numerical values of the energy supplied to the calo- 
rimeter by these transients are based on separate 
determinations of the transient voltage-time curves, 
as observed on an oscilloscope, for the power supplies 
operated in the appropriate mode of operation. Cor- 
rections to account for the slightly different average 
values of power supplied to the heater in the actual 
electrical calibrations were made assuming the length 
of time required to switch the power supply from the 
dummy to the calorimeter heater was constant. 
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Table 5. Summary of calculation of (— AE B ) 



Expt. 

No. 


Number 
expts. 


of 


Heater power 


Average £(cal),y°C" 1 


Standard deviation of an 
experiment,./ °C ' 


Standard deviation of 
average, J °C ' 


95% confidence limits 

y°c-« 


7-15 


9 

2 
4 

1 


110W (constant voltage 
60 W (constant voltage).. . 
270 W (constant current I. 
300 W (constant voltage).. 


14,127.79 
14,127.72 
14,125.19 
14,128.36 
14,127.79 


0.99 (0.007%) 

0.99 

1.6 


0.33 
0.73 
0.80 

0.30 


0.76 


16-17 


1.60 


18-21.... 


2.55 


22 




Adopted value.. 


0.99 (0.007%) 


0.68 









Summary of benzoic acid combustions 



Expt. No 



Number of 
expts. 



Average 



E(cal)/(-A£„),g°C- 



Standard deviation of an 
experiment, g °C~' 



Standard deviation of 
average, g °C " ' 



■> confidence limits 

g c-> 



(0.0034%) 



Calculation of— AE B 

Nominal (— A£«) 

Correction for R aMi 

Correction for F 

Corrected (— LE H ) — 

with an uncertainty interval (random) of 



26,433.69 J g 
+ 0.21 j g" 
+ .15 J g- 
26,434.05 J g 
17 Jg-> 



0.0000193 



5.4. Summary of results 

The calculations of the energy equivalent of the 
calorimeter are summarized at the top of table 5. 
Ninety-five percent confidence limits are calculated 
by multiplying the standard deviation of the mean by 
the appropriate factor of the Student ^-distribution. 
Assuming equal precision, standard deviations com- 
puted for experiments 7-15 and for 16-17 are pooled 
together. Comparison of the average values of ^(cal) 
listed in table 5 for the three heater powers used in 
experiments 7 through 21 show they lie within their 
combined uncertainty (random) intervals. However, 
the results of experiments 18-22 were not used in 
computing a weighted average for ^(cal) for the fol- 
lowing two reasons. 

First, it will be noted from columns 2 and 4 of table 
4 that the contribution of Qr to /^(cal) is 0.2 percent in 
the case of experiments 18-21. For the other experi- 
ments, the power supplies were operated in the con- 
stant voltage mode and the contribution of Qr to £(cal) 
is less than 0.0012 percent. Since we feel that a sys- 
tematic error in the estimate of Qr of 20 percent is not 
unreasonable, the difference between the average 
^(cal) of experiments 18-21 and the remainder is 
probably caused by a systematic error due at least in 
part to the error in Qr. 

Second, it will be noted that experiments 18-22 
were carried out with the highest heater power and, 
thus, the shortest time necessary to warm the cal- 
orimeter 3 °C. This resulted in the unfortunate situa- 
tion that we *were able to make only a few measure- 
ments of heater current or voltage when the change in 
voltage, or current, was the most rapid. We feel the 
uncertainty introduced in ^(cal) by the uncertainty 
in our estimates of the time dependence of the voltage, 
or current may be appreciable. Thus, the agreement 



of experiment 22 with 7-17 is regarded as being some- 
what fortuitous. Experiment 22 was given zero weight 
on the basis that we suspect its uncertainty is con- 
siderable larger than experiments 7-17. 

The results of the benzoic acid combustions and the 
value determined for AE B for benzoic acid are sum- 
marized at the bottom of the table 5. The precision 
of kE B was calculated as suggested by Rossini [27]; 
the contribution of the imprecision due to the various 
factors taken as constants in determining E (cal) was 
neglected since it was estimated to be 10 ppm. 

. 6. Appraisal of Experimental Results 

6. 1 . Random Errors 

The standard deviation of a single combustion ex- 
periment with this particular calorimeter had never 
before been lower than 0.012 percent, for an extended 
series of measurements. It is clear from table 5 that 
the precision of the present measurements is sub- 
stantially better than the previous work. 

We were particularly interested to see whether or 
not we could account for the standard deviation of a 
measurement of either —E( cal)/ &E H or £(cal) in terms 
of the various factors we might expect to vary from 
experiment to experiment. 

Common to both sets of measurements were the 
uncertainty in the corrected temperature rise and the 
absence of a buoyancy correction to the weight of the 
water placed in the calorimeter can. 

The uncertainty in the corrected temperature rise 
was evaluated from the variances of the estimates of 
the various parameters appearing in eqs (2) and (3), 
using the usual propagation of error formulas [27] 
and assuming that the "average" root mean square 
deviation of 6(t) for the initial and final drift periods 
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may be treated as a random error. These "averages" 
were arbitrarily multiplied by V2 since they are prob- 
ably too small due to the use of eq (1) (see sec. 5.1a). 
The error in the integral in eq (5) was based on the 
approximations that the main error is contributed by 
the errors in the times recorded during the rapid tem- 
perature rise of the main period and that the integral 
may be computed by the trapezoidal rule. 

The uncertainty due to the absence of the buoy- 
ancy correction in the calorimeter weight was esti- 
mated from the known variation in air density during 
the benzoic and experiments (12 mg/1 7 ) and the ap- 
proximate volume of the calorimeter can plus water 
(2.59 1). The resulting percentage error was arbi- 
trarily multiplied by two because it is quite evident 
from the method used to adjust the weight of water in 
the can (see [3]) that the aforementioned variation 
in air density is a minimum estimate. Other factors, 
such as convection currents, may also make a sub- 
stantial contribution to the weighing error. 

The random error in £(cal)/A£« due to the uncer- 
tainty of the mass of the benzoic acid sample is based 
on a conservative estimate of ±10 /xg for a sample 
weight of 1.5 g. 

The estimates of the varous errors contributed by the 
parameters which make up Fi are based on the devia- 
tions of their measurements. The standard current 
resistor, Ro.oi, was assumed to have a stability of 2 ppm. 

The estimate of the random error in E(ca\) due to 
the random error in the integrals in eq (9) was based 
on the approximation that 13 smoothed pairs of values 
of Ev and Ej separated by 30 s time intervals were 
numerically integrated by the trapezoidal rule. The 
errors in the smoothed values of Ev and Ei were taken 
as ± 2 X 10~ 7 potentiometer units. 

The estimate of the random error due to the heater 
leads is based on the analysis given by Ginnings and 
West [28]. The dimensionless parameters L 8 and L c 
of that analysis were calculated to be 0.34 and 0.28, 
respectively. The variation in the heater lead resist- 
ance was estimated to be ±0.002 (1; the heater resist- 
ance is 32 (I. 

A summary of the errors and a comparison with the 
observed standard deviation of an experiment is given 
in table 6. The discrepancy between the observed and 



estimated values for the electrical calibrations is 
probably greater than the uncertainty of our estimates. 

6.2. Systematic Errors 

One of the main assumptions customarily made in 
using a stirred-water "isoperibol" calorimeter is that 
the correction to the temperature rise due to heat 
transfer from the calorimeter jacket and heat delivered 
by the stirrer can be computed from Newton's cooling 
law and temperatures measured inside the calorimeter. 

It is generally understood that this assumption would 
be valid if the calorimeter temperature were strictly 
uniform throughout the calorimeter. 8 Real calorime- 
ters have, of course, temperature gradients or, equiv- 
alently, temperature differences due to the effects of 
lags. Thus, the surface temperature of the calorimeter, 
which determines the rate of heat transfer between 
the calorimeter and jacket, is in general different from 
that measured by the thermometer, even though it is 
placed as close to the calorimeter wall as is practical. 

The effect of a temperature gradient on calorimetric 
measurements has been analyzed to some extent for 
"isoperibol" calorimeters [29, 30, 31], and in detail 
for aneroid adiabatic calorimeters [32]. Harper [29] 
has shown that if the lag of the temperature measure- 
ment system (i.e., including both thermometer and 
electrical detecting system) is the same throughout 
a calorimeter experiment then the lag causes no error 
in the corrected temperature rise. White [30], in 
effect, considers the temperature gradient to be the 
superposition of those due to (1) heat transfer from the 
jacket and stirring, and (2) heat released in the calo- 
rimeter during the main period. By an argument 
identical to that used by Harper, he showed that the 
lag due to the first of these sources causes no error 
in the corrected temperature rise, while the lag due 
to the second source will cause an error. However, if 
the latter lag is a constant which is independent of 
the source of heat (i.e., electrical heater or combustion 
bomb), then it will cause no error in the heat of com- 
bustion of benzoic acid. White felt that the constancy 
of the lag due to different heat sources could be veri- 
fied by appropriate experiments. 



7 Square root of the sum of squares of the deviations divided by the number of measure- 
ments minus one. 



8 Convective heat transfer between the calorimeter and jacket must also be negligible. 
For this discussion, we assume this is so. 



TABLE 6. Estimated random errors in — E(ca/)/( — AE B ) and E(cal) 



Benzoic acid 

(-£(cal)/AA'„) 



Electrical calibrations 
(-£(cal)) 



Corrected temperature rise 

Bouyancy correction 

Mass of benzoic acid 

[ EyE,dt 

klE' m 

£*, 

Ko.e, 

Volt box ratio! 

Heater leads 

Estimated standard deviation of a measurement 
Observed standard deviation of a measurement. 



0.0022 
0.0018 
0.0007 



0.0030 
0.0034 



0.0019 
0.0018 



0.0003 

0.0001 
0.0001 
0.0002 
0.0001 
0.0005 
0.0027 
0.0070 
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Table 7. Summary of values of (— AE B ) for benzoic acid 



Author 



Sample 



Value Uncertainty 



Jessup and Green [2. 18, 33] (1934) 
Prosen and Rossini [35] (1939-41).. 

Jessup [2.33] (1942) 

Coops et al. [31] (1946-8) 

Coops et al. [31] (1954) 

Challoneretal. [34] (1954) 

This work (1966) 



39d, 39e 

39e 

39e,39f, B 

39b, 39f, VUS, 

39b, VUS 

C 

39i 



Jg- 1 

26,432.0 
26,434.7 
26,433.8 
26,438.0 
26,435.0 
26,436.0 
26,434.0 



J g-> 



±2.6 

±2.2 
±2.6 

±4 
±4 
±4 
±3.3 



Coops et al. [31], have suggested one method of 
making the lag due to different heat sources the same 
and have incorporated this in the design of their 
calorimeter. 

An experimental test that has frequently been ap- 
plied is that the lag is a constant if it can be shown that 
the energy equivalent of the calorimeter is invariant 
with the power supplied to the electrical heater. This 
criterion has been shown to be incorrect by West [32] 
for a calorimeter in which heat flow equations are 
linear. It is interesting to note that part of this argu- 
ment depends upon showing that the contributions to 
the correction on the observed temperature rise of 
the temperature transients due to turning the elec- 
trical heater on and off exactly cancel. Indirect evi- 
dence that this cancellation is the case for our calorim- 
eter is given in section 7. This conformity of our 
calorimeter with West's analysis of behavior of tran- 
sients suggests that the remainder of West's argu- 
ment may be applicable to the Dickinson calorimeter. 
The invariance of the energy equivalent found in the 
electrical experiments with varying electrical power 
(see table 5), thus, does not provide information about 
possible systematic errors. 

A series of tests was carried out to see whether or 
not measurements based on the average surface tem- 
perature of the calorimeter in place of that registered 
by the platinum thermometer would give a different 
value of AEk for benzoic acid. On the basis of these 
tests, which must be considered to be preliminary in 
nature, it was found that there is a difference between 
the temperatures of the calorimater surface and the 
water in the vicinity of the resistance thermometer 
and that this difference varies with time in the elec- 
trical calibrations in a different way than in the benzoic 
acid experiments. Because of certain experimental 
difficulties, we assigned an uncertainty to the differ- 
ence in values of &E B calculated from these observa- 
tions that was about the same as the magnitude of the 
difference. Until a more reliable estimate based on 
more accurate tests is obtained, we have somewhat 
arbitrarily assigned an uncertainty of ±0.01 percent 
of AEb to take into account the possible correction for 
the effect of surface temperature. 

The net contribution of other sources of systematic 
errors in the electrical calibrations is estimated to be 
at most ±0.004 percent of hE B . 



6.3. Discussion of Results 

In table 7, we have listed the more recent determi- 
nations of (—AEb) for benzoic acid in joules per gram 
of sample. 9 Samples indicated by the number 39 fol- 
lowed by a letter refer to batches of the benzoic acid 
standard sample issued by NBS. B refers to a sample 
of 39e purified by fractional crystallization from 
benzene. VUS and C refer to benzoic acid samples 
from sources other than NBS and are described in the 
corresponding references. 

Uncertainty intervals listed in table 7 are those cited 
by the authors except for the work of Jessup and 
Green [2, 18, 33] which we assumed to be equal to 
the later work of Jessup [2] . The bounds for the un- 
certainty interval for our own work was computed by 
the procedure recommended by Rossini [27] as 



V(1.7) 2 +(l.l) 2 +(2.6) 2 . 

Aside from our own work, only that of Challoner, 
Gundry, and Meetham [34] incorporates an uncer- 
tainty due to the effect of surface temperature. In the 
latter, a correction of — 0.046 ±0.0057 percent had to 
be applied to the energy equivalent of the calorimeter 
to correct for the difference in surface temperature 
effects during the main period of an electrical cali- 
bration and a benzoic acid combustion experiment. 
Our preliminary experiments suggest that the magni- 
tude of the surface temperature correction is not 
insignificant for our results although it is smaller than 
that of Challoner et al. A substantial contribution to 
the uncertainty of our results has been allowed be- 
cause of the unknown magnitude of our surface tem- 
perature correction. It would be preferable, in spite 
of the difficulty of measuring it, to evaluate the cor- 
rection experimentally in absolute calorimetric meas- 
urements unless the calorimeter has been specifically 
designed to eliminate the error. 



9 We bave been informed by private communication from Mosselman and Dekker [36], 
and from Head [37], of very recent electrical determinations of the energy of combustion 
of benzoic acid in their respective laboratories. These determinations are both based on 
measurements made in calorimeters containing water flow channels of the type devised by 
Coops et al. [31], and both are reported to be in excellent agreement with the value we list 
for this work in table 7. The work of Mosselman and Dekker includes some measurements 
on NBS Standard Sample 39i. 
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7. Appendix. Analysis of Main Period of the 
Electrical Calibration Experiments 

If the temperature transients due to turning the elec- 
trical heater on and off make no net contribution to 
the corrected temperature rise, the two areas between 
the main period curves of the electrical calibrations 
and the straight lines given by eqs (10a), (10b), and 
(10c) should be equal since the heater power does not 
vary with time. 

0(t) = 6 i +(d0/dt) i(t-tj); ti^t^h (10a) 

0(0 = 0:*+ {d6/dt) 3 (t — t 3 ); t\ ^t^H + tb (10b) 

0(0 = 0/+ (dOldt)f(tj-t) ; t 2 + t b ^ t =s= t f . (10c) 

Equations (10a) and (10c) are linear extrapolations of 
the initial and final drift period curves into the main 

period, and the constants 6 :i , tz and \~j~) of eq (10b) 

are determined from a linear least squares fit of the 
main period data during the rapid temperature rise. 
The times t\ and £ 2 + tt> are the times of intersections 
of eq (10b) with eq (10a) and (10c). 

To test whether or not the areas are equal, we noted 
that if they are, then the time t\ should be approxi- 
mately equal to time 10 fa and both should be inde- 
pendent of the power supplied to the heater or (dO/dt)*. 
Comparison of columns 3 through 5 of table 8 shows 
this is very nearly true. Further, the equality of the 
areas means that the integral in eq (5) for t x is the area 
enclosed by the lines of eqs (10a) through (10c) and 0,. 
Corrected temperature rises computed using values 
tf—tjc determined in this manner differed by less than 
one part in 10 5 from those listed in section 5. The 
latter are based on values of t/—t x determined by 
the more accurate spline integration procedure. 



Table 8. Average main period characteristics, electrical calibration 



Heater power 
watts 


No. of 
Expts. 


min 


min 


(dd/dt)* 
°C min -1 


min 


'max 

min 


[28] 
[29] 
[301 


60 
110 
270 


2 
9 
5 


0.238 
0.242 
0.230 


0.224 
0.235 
0.254 


0.252 
0.480 
1.15 


11.2 
5.8 
2.5 


12.5-13 

8.5 
4-5.5 
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